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Abstract SVOM (Space-based multi-band astronomical Variable Object Mon- 
itor) is a future Chinese- French satellite mission which is dedicated to Gamma- 
Ray Burst (GRB) studies. Its anti-solar pointing strategy makes the Earth 
cross the field of view of its payload every orbit. In this paper, we present the 
variations of the gamma-ray background of the two high energy instruments 
aboard SVOM, the Gamma-Ray Monitor (GRM) and ECLAIRs, as a function 
of the Earth position. We conclude with an estimate of the Earth influence on 
their sensitivity and their GRB detection capability. 

Keywords Gamma-ray background • Earth occultation • Monte-Carlo 
O simulation ■ Gamma-Ray Burst 
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1 Introduction 

" After forty years of studies, we are still far from a full understanding of 

Gamma- Ray Bursts (GRBs). Accurate measurements of GRB parameters, 
such as their position, redshift, peak-energy and so on, are needed to fur- 

. thcr understand GRBs themselves and their use as astrophysical tools. SVOM 
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(Space-based multi-band astronomical Variable Object Monitor) [I] [5] is a fu- 
ture Chinese-French satellite mission which is dedicated to GRB studies. The 
payload aboard SVOM consists of two high-energy instruments, the GRM and 
ECLAIRs, together with two low-energy instruments. GRM is a Gamma Ray 
Monitor, and ECLAIRs is a coded-mask camera for X- and Gamma-rays. To- 
gether, they can provide GRB localizations and GRB spectral observations 
from a few keV to a few MeV. This wide spectral range is crucial to determine 
as accurately as possible the GRB peak-energy, which is a key parameter to 
make use of GRBs as standard candles and extend the measurement of the 
cosmological parameters at large redshifts pQ [3] ■ 

SVOM is a Low Earth Orbit (LEO) mission with an altitude of 600 km 
and an inclination of 30°. Its instruments point close the anti-solar direction 
during a large fraction of the orbit in order to permit the follow-up observations 
by large ground-based telescopes to optimize the redshift measurements [I] [4] . 
This strategy makes the Earth get in and out of the detectors Field of the View 
(thereafter FoV) once per orbit. The detectors will not operate when crossing 
the South Atlantic Anomaly (SAA) for protection from the too high particle 
flux, which results in a decrease of the observation time. These factors will 
cause the instrumental background to change when the Earth gets in and out 
of the FoV and with time. In this paper, we present a study on how the gamma- 
ray background of GRM and ECLAIRs changes with the Earth position in the 
FoV. The impact of these changes on the sensitivity and GRB detection rate 
are estimated as well. This background research laid the groundwork for the 
investigation of the GRB trigger strategy of GRM and ECLAIRs. This paper 
is organized as follows: 

In this section, we briefly describe the two high energy instruments 
aboard SVOM. 

USlThe detector mass-models, the simulation method and the spectral mod- 
els of gamma-ray background sources are described in this section. 

<2] The background simulation results, the computations of the detector 
sensitivities and the GRB detection rates are presented in this section. 

<J3] We conclude with a concise summary and a discussion of some limita- 
tions of our method as well as some further improvements. 



2 Instruments 

The main parts of GRM are two identical Gamma Ray Detector (GRD) units, 
which are spectrometers without localization capabilities, sensitive in the en- 
ergy range from 30 keV to 5 MeV. Each GRD is equipped with one auxil- 
iary Calibration Detector, which mainly consists of a plastic scintillator and 
an embedded Am source, in the collimator above the entrance window. A 
Charged Particle Detector is mounted to alert the high-voltage supplies of 
PMTs to switch off or on. The GRD consists of three scintillator layers (see 
Figure[T]) which are plastic scintillator, Nal(Tl) and CsI(Na) respectively. They 
are glued together and viewed with the same light guide coupled to a PMT. 
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The pulse shapes generated by the different scintillators can be distinguished 
by the Pulse Shape Discrimination circuit (PSD). The 15 mm thick Nal(Tl) 
crystal with a diameter of 190 mm works as the main detecting element. The 
CsI(Na) crystal, which is 35 mm thick and has the same diameter as Nal(Tl), 
located behind Nal(Tl) is another important detection element and serves also 
as ant i- coincidence element against photons from behind. The element located 
on top of Nal(Tl) is the organic scintillator, dedicated to reject background 
events due to charged particles. The structure of organic scintillator has been 
modified since the previous design described in [5] cannot shield-off efficiently 
the photons from the side and its design is more complex. The beryllium plate 
with a thickness of 1.5 mm is chosen as entrance window of this triple phoswich 
detector. A collimator made of tantalum is located in front of the scintillator 
case to reduce the background by limiting the FoV to 2.5 steradians. The main 
components of GRD are shown in Figure [TJ 



Collimator 
(Tantalum) 



Nal, 15mm 




1.5mm thick Be 
k 6mm thick PS 



Be,PS,Nal,Csl 



Diameter: 190mm 



Bracket 



Front-End Electronics 



Fig. 1 View of the main components of GRD. 



ECLAIRs is a coded-mask imaging camera with a 2 sr-wide FoV, with a 
source localization accuracy better than 10 arcmin, sensitive in the 4 keV - 250 
keV energy range. Its main entity is an imaging Camera for X- and Gamma- 
rays (CXG) which mainly consists of a detection plane of 80x80 CdTe semi- 
conductor detectors, a coded-mask located above the detection plane, a multi- 
layer lateral shield and mechanical structures, etc. ECLAIRs is also equipped 
with an onboard Scientific Trigger Unit (UTS) [B] which seeks in real-time 
for GRBs and determines their localization for follow-up observations. For a 
detailed description of ECLAIRs, see [7] [5]. 
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3 Monte-Carlo simulations of the instruments gamma-ray 
background 

In this section, the mass models of GRM and ECLAIRs used to perform our 
simulations, the simulation method as well as the corresponding background 
models are described. 

3.1 Mass models and simulation method 

The simulations have been performed with release 4.9.3p02 of the Geant4 
toolkit [9 . In accordance with the design of GRM, whose components and 
main characteristics are described in Section [2J its mass model was built. In 
order to be close to reality, the other payloads, the pay load board and the 
satellite platform are also considered in the mass model, as shown in Figure 
IH(left). But currently only GRDs are detailed, the others are represented by 
shells of aluminum with shapes and thickness corresponding to their latest 
design. 




Fig. 2 Schemes of the mass models of GRM and ECLAIRs. (left) Scheme of the mass 
model of GRM (see Section[2]and Figure^ as well as the shells of other payloads. The grey 
boxes correspond to the satellite platform and payload board (not the exact structure of 
the spacecraft), (right) Scheme of the mass model of CXG as well as the cold plate and the 
readout electronics box placed under the camera. The gray box corresponds to the satellite 
platform. 

For ECLAIRs, the CdTe detectors, the ceramics, the ASICs, the multi- 
layered shield composed of Pb/Cu/Honeycomb/Cu, the coded mask made of 
a Ta-W alloy, the multi-layer thermal coating insulation (MLI) above the coded 
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mask, the TA6V upper and lower supports of the mask, as well as simplified 
versions of the cold plate in AlBeMet and the readout electronics boxes are 
considered in the mass model. The platform is represented by a 1 m 3 and 300 
kg pseudo-aluminum cube. 

In the background simulations, the most important point is to study how 
the background varies with the Earth position in the FoV. In practice, we keep 
the satellite static and let the Earth rotate around it. The establishment of 
the coordinate system is shown in Figure [21 the Z axis points to the opposite 
direction of the FoV and the coordinate origin is at the center of the satellite. 
When the Earth rotates from Z to -Z, the angle changes from 0° to 180°. Wc 
divide the half orbit equally into 7 pieces, with an angle-interval of 25.7°. For 
each piece of orbit, the Earth position can be represented by a typical position 
approximately. And the zenith angle of each typical position in the coordinate 
is used to express the Earth position. The typical Earth position of Rotate 0° 
means that the Earth is on Z and Rotate 180° is the Earth position which is on 
-Z. Other typical positions in between are Rotate 30°, Rotate 60°, Rotate 90°, 
Rotate 120° and Rotate 150° respectively. When the Earth rotates from -Z to 
Z, this half orbit is also divided equally into 7 pieces as the same as that from 
Z to -Z. Hence the 7 pieces from -Z to Z are axisymmetric with those from Z to 
-Z respectively, and can still be represented by the 7 typical positions (Rotate 
0° - Rotate 180°). Here we assume that for different Earth positions with the 
same zenith angle but different azimuth angles, the background level is the 
same. With the assumptions and simplifications above, we can now simulate 
the background levels with the 7 typical Earth positions, and each position 
lasts 1/7 of the orbital period. These assumptions will be used in Section FOl 
to derive a first order estimation of GRB detection rate. 



3.2 Space environment 

The satellite will be subject to different sources of background in space: (i) 
cxtragalactic components including cosmic X- and gamma-rays, protons, elec- 
trons, positrons, and cosmic He-nuclei, etc.; (ii) near Earth components con- 
sisting of albedo gamma-rays, albedo neutrons, albedo electrons, albedo positrons 
and the secondary protons located under the radiation belt; (iii) as well as the 
particles from the Sun. 

In this paper three kinds of gamma-ray background are investigated in de- 
tail, which are the main background components in the enregy range of our 
interest and include cosmic gamma-rays, albedo gamma-rays and reflected cos- 
mic gamma-rays. We adopt the latest spectral models which are Moretti2009 
|10) . Sazonov2007 [11 and Churazov2006 [I2][l3] respectively. For convenience, 
we will denote cosmic gamma-ray background, reflected cosmic gamma-ray 
background and albedo gamma-ray background as CXB, REFLECTION and 
ALBEDO respectively in the following context, and if there is no specific dec- 
laration the term background only means the gamma-ray background. 
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For ECLAIRs, the background is dominated by the gamma-ray background 
below 100 keV outside SAA. And the proportion of the background due to the 
direct interactions with particles such as electrons, protons and neutrons is 
about 10% [7] [T3]. Concerning the delayed background mainly introduced by 
the trapped particles inside SAA, it is outside the scope of this work and we will 
not discuss it in detail. For GRM, in addtion to our simulations of gamma-ray 
background, G. Li estimated the background induced by the cosmic protons 
and electrons, the secondary protons and electrons, as well as the trapped 
protons in SAA, when the Earth is at Rotate 0° using the method provided in 
[IS] . With these two parts of work, we obtained that (G. Li, personal commu- 
nications, 2011) the average proportion of the gamma-ray background in all 
kinds of background is approximately 80% in the energy range 50 keV-300 keV 
for GRM_NaI. The average proportion of the delayed background caused by 
the trapped protons in SAA region is less than 15% below 300 keV, and is ap- 
proximately 50% above 550 keV. For GRM_CsI, the proportion of the delayed 
background is slightly higher. The background which is used to do calculations 
in Section [4.31 only includes the three kinds of gamma-ray background. 

3.2.1 Cosmic Gamma-ray Background model 

The spectral model of Moretti2009 which is described by Equation (1) is 
adopted for CXB simulations. This model was produced by taking advantage 
of both Swift-XRT and Swift-BAT databases in the energy band 1.5-200 keV. 
The spatial distribution is isotropic in a 4tt solid angle. The Earth occultation 
is considered in the post-processing after the simulations with Geant4. 

(i) 

arj \ 29.0 J +^ 29.0 - 1 

Ajello2008 [TB] is another recent model with a smoothly joined double 
power law. It was obtained by fitting the data in the 2 keV-2 MeV range. 
Below 500 keV, the relative error is less than 20% using model Moretti2009 
and Ajello2008. And above 500 keV, the statistical fluctuation is much larger 
than the relative error. Other models such as Sreekumarl998 ]T7] and Gru- 
berl999 [TH] were also compared. Since there is not much difference using 
these models for GRM, and considering that Moretti2009 is a credible model 
for ECLAIRs, we finally determined to use the model Moretti2009 for both 
GRM and ECLAIRs. 

3.2.2 Reflected cosmic gamma-ray model 

The outer layers of the Earth's atmosphere reflect part of the incident X- and 
gamma-ray spectrum via Compton scattering. The spectrum in the 1-1000 keV 
energy range reflected by the Earth's atmosphere was calculated by Churazov 
et al. (2006). The ratio of the reflection spectrum to the incident spectrum 
is approximated by Equation (2). The incident spectrum is the CXB model 
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suggested by Gruber et al. (1999). The reflection is most significant near 60 keV 
and declines towards lower or higher energies. Since the incident directions of 
photons are connected with Earth positions, some methods used in simulations 
are similar to those for Albedo photons described in Section f3. 2. 31 



A(E) = 



1.22 



2 -93 +(^8) 4 v 0.123 +(^3) 



E \3.44 



MM 

\28.5> 



-2.54 



+ (HO) 1 ' 57 - 0-37 l + (^Y 



1 + (-2-) 

± ^ V 91.83/ 



3.44 



(2) 



3.2.3 Albedo gamma-ray model 



The albedo gamma-ray background consists of the atmospheric hard X-ray 
emission which is induced by cosmic rays such as energetic protons, alpha 
particles, cosmic electrons and positrons as well as secondary electrons and 
positrons. 

Sazonov et al. found that the shape of the spectrum emergent from the 
atmosphere in the range 25-300 keV is mainly determined by Compton scat- 
tering and photo-absorption, and is almost insensitive to the incident cosmic- 
ray spectrum. They provided a fitting formula for the hard X-ray surface 
brightness of the atmosphere which would be measured by a satellite-borne 
instrument, as a function of energy, solar modulation level |19j . geomagnetic 
cut-off rigidity R cut and zenith angle 0; see Equations ([3]), dU) and (|5|). 

= ( I \-^ C ( E [ph/cm 2 /s/sr/keV] (3) 

\4iheV) ~ t[ .44keV> 



C = 3M1 + M) 1-47 x 0.0178 x [(^2.8^ + ^ /eTO a /a/ar] 



[ 13< | > U.25( 1 + 2.5 < |,0.4)J 



(4) 



where, fi = cost 



Rcut = 59 ' 4 X [l + (l + c OS 3A m 2^<)VT [Gy] (5) 

Equation allows us to predict the hard X-ray surface brightness of the 
atmosphere or its part when observed from a spacecraft, provided that its 
orbit is higher than 50 to 100 km. As a practical example, Sazonov et al. used 
this formula to predict the hard X-ray flux which would be measured by a 
satellite-borne instrument whose FoV covers the entire terrestrial disc. In our 
case, SVOM is a LEO mission with an altitude of 600 km and an inclination 
of 30°. If we suppose that SVOM is aiming towards the nadir direction, only 
a part of the atmosphere is visible by the instruments. Taking into account 
that mainly the upper part of the atmosphere is producing the hard X-ray 
emission, we consider that the albedo source is at an altitude of 100 km. With 
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a satellite altitude of 600 km the projected radius of the visible atmosphere 
spherical cap is ~2400 km, see the region circled by white line in Figure [31 

For a given position of the satellite on the orbit, we have to locate the 
atmosphere spherical cap visible to the instruments. Inside this spherical cap, 
the magnetic latitude A m varies and therefore the geomagnetic cut-off rigidity 
is not uniform, as shown in Figure [3J For different satellite positions, we do 
not intercept the same area of the geomagnetic rigidity map. As a result, the 
global ALBEDO level varies with the satellite position along the orbit. 




Fig. 3 Geomagnetic cut-off rigidity (unit: GV) map. The region surrounded by a white line 
is the atmosphere spherical cap visible by SVOM when the satellite position is 30° N, 287° E. 

Hard X-ray surface brightness use 

In order to compute the hard X-ray surface brightness at a given satellite 
position we propose the following methodology: 

For a given satellite position, we define a new reference frame OXYZ in a 
spherical coordinate system with the Z axis going from the Earth center trough 
the SVOM satellite. The atmosphere elementary elements of the spherical 
cap are represented by the coordinates 9\ and (pi in the OXYZ frame (see 
Figure 4.). In practice we divided the atmosphere spherical cap into 2400 
elementary areas with 20 9\ and 120 ifi steps (a compromise between accuracy 
and computation time). For each elementary area, we do the calculations as 
follows: 

- We determine the geomagnetic latitude A m by performing a reference frame 
change from the OXYZ frame to the geomagnetic reference frame. The 
reference year for the geomagnetic reference frame is 2015 when the Earth's 
magnetic South Pole is located at (80.367°N, 72.624°wj±]. 



1 http:/ /wdc. kugi.kyoto-u.ac.jp/igrf/gggm/ 



Title Suppressed Due to Excessive Length 



9 



- Knowing A m , we calculate R cut which depends on the zenith angle 9 and 
on the azimuth angle £, see Equation ([3]). Following the simplification pro- 
posed by Sazonov et al. we assume £ = 0. 

- We compute the brightness according to Equation Qj. And the surface- 
brightness map corresponding to the satellite position (30°N, 287°E) is 
represented in Figure [SJ 



Fig. 4 Schematic drawing of albedo gamma-ray background source model, (left) Three- 
dimensional schematic drawing of the albedo gamma-ray background source model. The 
fan-shaped area is the elementary element defined by 9\ and ip±. The position of (6i, 
</?i=0) corresponds to the point of (satellite latitude - 0\, satellite longitude), (right) Plane 
schematic drawing of the albedo gamma-ray background source model (the straight line 
around SVOM represents the radiation source in the simulation model built in Geant4. In 
this work, the line P-SVOM is the borderline of the photon directions). 

For a given satellite position, by summing the contribution of each cell, we 
can compute the global hard X-ray surface brightness (ph/s). The highest level 
is obtained at the satellite position (30°N, 287°E) and the lowest level is ob- 
tained at (7°N, 151°E), see Figure [5] (left). Moreover, the brightness decreases 
as the solar modulation <P increases from 0.25 GV to 1.5 GV, see Figure [6] 
(right). In order to study the variation of the ALBEDO as a function of the 
Earth position, we decided to choose the most unfavorable case: we fixed the 
satellite position at (30°N, 287°E) with solar modulation 0.25 GV. 

Albedo gamma-rays computation 

For the chosen satellite orbit position, the above computations enabled us 
to find out the value of all parameters in Equations ([3][5]) and to compute 
a surface-brightness map providing us, for each elementary atmospheric cell 
of center PI characterized by coordinates (9i,<pi) and area A(9i,ipi), with a 
photon spectrum jrjsffiutpij E) in unit of photons/cm 2 /s/sr/keV. 
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Fig. 5 Hard X-ray surface brightness (unit: 10 ^ph/cm 2 /s/sr) at the satellite position 
(30°N, 287°E). 




SolsrMwJuiationlGV) 



Fig. 6 (left) The distribution of the global hard X-ray surface brightness along the satellite 
orbit. The highest level position is (30°N, 287°E), and the lowest level position is (7°N, 
151°E). (right) Variations of the global brightness as a function of the solar modulation 
when the satellite is at (30° N, 287° E). 



Considering that our Geant4 model of the SVOM satellite is enclosed in a 
sphere, we can compute the half-angle a of the cone subtended by the satellite 
from the center of each cell (see Figure HI right), and thus its solid angle 
fl(Q\,tp\) — 27r(l — cosa). 

After choosing an exposure time, with the spectrum 4§ (9\ , ifii ; E) , the solid 
angle f2(6i,ip\) and the area A(#i,</?i) of each elementary atmospheric cell, 
we can compute the total number of photons to shoot and the relative weight 
of each cell. Finally, we can generate the appropriate number of photons for 
each cell, randomly giving them the appropriate energy, direction and position. 
Given the chosen satellite attitude - the position of the Earth in the FoV, we 
apply a suitable coordinate transformation before injecting these photons into 
our Monte-Carlo model. 
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4 Results from gamma-ray background simulations 



Considering the sensitive energy ranges of instruments, the input energy spec- 
tra are extended to [10 keV, 100 MeV] for GRM and to [1 kcV, 100 MeV] for 
ECLAIRs in the simulations. The input spectra of CXB and REFLECTION 
were obtained directly from their energy spectral equations. For ALBEDO, 
the spectrum in Figure [7] is the average spectrum from the whole atmosphere 
spherical cap visible by SVOM when the satellite position is (30°N, 287°E) 
and solar modulation is 0.25 GV. Comparing the three spectra, we can see 
that the CXB flux is the highest at low energy (<~100 keV), the ALBEDO 
flux becomes highest as energy increases, and the REFLECTION has an in- 
termediate level at low energy. These features imply that the ALBEDO will 
play a very important role in the high-energy range. The difference between 
these three kinds of background is not only the generation mechanisms but 
also the spatial distributions. This suggests their influence on the final in- 
strument background will be varying with the Earth position. In this section, 
the instrumental backgrounds of GRM and ECLAIRs induced by these three 
kinds of background sources are described. In Figure [9j [TT] and Figure [13l we 
only show the results corresponding to Rotate 0°, Rotate 90°, Rotate 120° and 
Rotate 180°, other Earth positions lead to very similar spectra. The counting 
rates of each background component in various energy ranges corresponding 
to the Earth positions of Rotate 0° and Rotate 180° are described in Tabled] 
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4.1 GRM gamma-ray background 

GRM has two detecting elements Nal(Tl) and CsI(Na). Their backgrounds are 
shown and discussed respectively. Figure |5J and Figure |H] present the changes 
of the three kinds of background in the Nal as well as the total background 
with Earth positions. 

From Figure |5J the CXB reaches the maximum level when the Earth is at 
Rotate 0° and reaches the minimum when the Earth is at Rotate 180°. On 
the contrary, the ALBEDO and REFLECTION reach their maximum level 
when the Earth is at Rotate 180°. Figure [J2 shows that the CXB decreases 
and the other two backgrounds increase quickly and obviously when the Earth 
rotates more than 90°. Finally when the Earth gets to Rotate 180°, the CXB 
becomes the lowest background component at a relatively low energy range 
(less than 200 keV), while the ALBEDO is always the highest at relatively 
high energies (more than 300 keV). That is because its flux is always the 
highest one at energies above 100 keV. The REFLECTION, which is always 
the lowest background component above 200 keV, starts to become higher 
than the ALBEDO below 50 keV when the Earth starts to enter the FoV. 

From Figure [51 for each separate background we can see obvious differences 
between Earth positions only at relatively low energies (below 300 keV) instead 
of the whole energy band. That is because the shields around the detector are 
not as effective for high-energy photons as for low-energy photons. For the 
same reason we cannot see obvious differences in the background when the 
Earth rotates less than 90° from the back of GRM. 

The emission peak of 57 keV in the CXB spectrum is due to the element 
Tantalum in the collimator. When the Earth enters the FoV, less cosmic pho- 
tons can enter Nal directly through the FoV, and relatively the proportion of 
counts induced by the photons, which are generated by Compton scattering 
of the relatively high energy photons with Csl or other surroundings, becomes 
larger at low energy. Overall, Earth occultation makes the counts at low en- 
ergy decrease, and the photons which interact with Tantalum in the collimator 
don't vary significantly. Consequently, the proportion of the counts caused by 
Tantalum fluorescence becomes higher and the peak becomes more obvious as 
the Earth enters the FoV. The 511 keV peak caused by electron-positron anni- 
hilation is only obvious in the ALBEDO spectra. That is because the photons 
with energy more than 1022 keV from the ALBEDO source are much more 
numerous than those from the other two background sources. The 511 keV 
peak is detectable in the total background spectra, and can be used for energy 
calibration. 

The total background is shown in Figure [5] (right bottom). At low energy 
(less than 50 keV) , the background level is highest when the Earth is at Rotate 
0°. At an energy above 50 keV, the highest background is obtained when the 
Earth is at Rotate 180°. 

Figure [TO] and Figure [IT] present how the background of Csl changes with 
Earth positions. In general, the background of Csl is higher when the back- 
ground source is outside the FoV than that when it is inside the FoV. The 
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ALBEDO and REFLECTION components reach their maximum (minimum) 
levels when the Earth is at Rotate 60° (Rotate 180°) in 60 keV-250 keV. And 
in this energy range, the CXB component reaches the maximum (minimum) 
level when the Earth is at Rotate 180° (Rotate 60°). That is because of the 
design of Csl and its position inside GRM (see Figure [TJ. On one hand, Nal 
which is in front of Csl can reduce the background of Csl from the FoV, and 
the structures such as the front-end electronics and so on which are behind Csl 
can decrease the background from behind. One the other hand, the interaction 
cross-section of Csl is the largest when the photons come from the direction 
of Rotate 60°. See Figure 11, ALBEDO and REFLECTION decrease as the 
Earth moves close to the FoV when the CXB increases gradually. Finally, the 
CXB becomes the highest background at relatively low energies (below ~200 
keV). The ALBEDO is always the highest background at relatively high en- 
ergies (above ~200 keV). The REFLECTION is the lowest component in the 
whole energy range. And the total background of Csl is shown in Figure ITD1 
(right bottom). 

4.2 ECLAIRs gamma-ray background 

The same background models are applied to ECLAIRs and the results are 
shown in Figure Q2] and Figure 1131 The results are analyzed in the same way 
as what has been done for GRM. 

From Figure [T^l the CXB reaches its maximum level when the Earth is 
at Rotate 0° and it reaches its minimum level when the satellite faces the 
Earth. The ALBEDO and REFLECTION components increase gradually as 
the Earth moves close to the FoV. Figure [13] shows the CXB becomes lowest 
at relatively low energies when the Earth reaches Rotate 180°. The ALBEDO 
is always the highest in the relatively high energy range. The REFLECTION 
plays a more and more important role when the Earth moves near the FoV. 

There are a few peaks in the background spectra due to element emission 
lines, including the peaks of 8 keV and 75 keV from the element Copper 
and Lead in the shield respectively, and the peak of 57 keV from the element 
Tantalum in the mask. The most intense Cu-Ta-Pb lines will be used for energy 
calibration. 

The total background is presented in Figure [T^] (right bottom) . The Earth 
occultation results in the background with marked difference at low energy 
(below ^50 keV). And the highest background level is obtained when the 
Earth is at Rotate 0°. At energies above 100 keV, the highest background 
level is obtained when the Earth is at Rotate 180°. 



4.3 Impact on the sensitivity and the GRB detection rate 

The background variations for GRM and ECLAIRs with the relative Earth 
positions defined in Section 13.11 are presented in Section 14.11 and Section 14.21 
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Fig. 8 The gamma-ray background levels of GRM_NaI change with the Earth positions, 
(left top) CXB of GRM_NaI changes with the Earth positions, (right top) REFLECTION of 
GRM_NaI changes with the Earth positions, (left bottom) ALBEDO of GRM_NaI changes 
with the Earth positions, (right bottom) The total gamma-ray background level of GRM_NaI 
changes with the Earth positions. In Figure [8] - Figure 1 131 the meanings of different lines 
shown in the legend on the corresponding top right panel are valid for all the four panels. 



Title Suppressed Due to Excessive Length 



15 



GRMJMal, Rotate 0*deg 

— i | r 



GRMJMal, Rotate 90*deg 

-TTTT| I T 




Energy(keV) 
GRM_Nal, Rotate 120*deg 



Energy(keV) 
GRM_Nal, Rotate 180*deg 




10' 



10 _J r 



> 

^io- 4 



O ln -5 



1CP f 



10' 



10 

























A., , ! 



10 10 
Energy(keV) 



10 10 
Energy(keV) 



Fig. 9 Three kinds of gamma-ray background of GRM_NaI change with the Earth positions. 



Based on these simulations, the peak fluxes at detection threshold in the energy 
range 1-1000 keV were computed using the method described in [20] . In the 
computation, we considered two BAND spectra [21] of GRBs: one spectrum 
with a = -0.5, j3 — -2, the other one with a = -1, /3 — -2. GRBs located on 
the axis of the FoV, and the background level with Earth position Rotate 0° 
was used. The integration time was 1 s and the significance a was 5.5. The 
4-50 keV and 50-300 keV energy bands are selected respectively for ECLAIRs 
and GRM considering their efficiency. The results are shown in Figure 1141 
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Fig. 10 The gamma-ray background levels of GRM_CsI change with the Earth positions, 
(left top) CXB of GRJVLCsI changes with the Earth positions, (right top) REFLECTION of 
GRM.CsI changes with the Earth positions, (left bottom) ALBEDO of GRM.CsI changes 
with the Earth positions, (right bottom) The total gamma-ray background level of GRM_CsI 
changes with the Earth positions. 



In general, ECLAIRs is sensitive to the GRBs with low Epeak and slightly 
more sensitive to the harder GRBs with a — -0.5 when Epeak < 100 keV. GRM 
is sensitive to GRBs with relatively high Epeak, especially to the harder bursts 
with a = -0.5. As shown in Figure PHI comparing with ECLAIRs, GRM_NaI 
is more sensitive for bursts with Epeak > 50 keV when a = -0.5, and more 
sensitive for bursts with Epeak > 90 keV when a = -1. GRM_CsI is sensitive 
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Fig. 11 Three kinds of gamma-ray background of GRM_CsI change with the Earth posi- 
tions. 



to the hard GRBs with high Epeak. The sensitivity of GRM and ECLAIRs are 
complementary which makes SVOM sensitive in a wide Epeak range, which 
allows accurate Epeak measurements. 

Based on the background-simulation results in the previous section, the 
corresponding sensitivities (Tabled and GRB detection rates fFigure [T5|) were 
calculated in the specific energy range with the given significance considering 
the GRB trigger efficiency, the observation-time loss due to SAA passages, etc. 
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Fig. 12 The gamma-ray background levels of ECLAIRs change with the Earth positions 
(left top) CXB of ECLAIRs changes with the Earth positions (right top) REFLECTION of 
ECLAIRs changes with the Earth positions (left bottom) ALBEDO of ECLAIRs changes 
with the Earth positions (right bottom) The total gamma-ray background level of ECLAIRs 
changes with the Earth positions. 



In Table [2j the FoV is the Field of View which is not obstructed by the 
Earth; the Sensitivity, which is the minimum flux that can be detected, is 
derived from the formula ([5]). For ECLAIRs, the Earth occultation results in 
the background with magnitude difference at low energy (see Figure I12|) . This 
induces the value of sensitivity (1.46 ph/cm 2 /s) corresponding to Rotate 120° 
less than that (1.62 ph/cm 2 /s) corresponding to Rotate 0° (see Tabled]). This 



Title Suppressed Due to Excessive Length 



19 




Fig. 13 Three kinds of gamma-ray background of ECLAIRs change with the Earth posi- 
tions. 



means that the appropriate Earth occultation can improve the sensitivity of 
ECLAIRs. Taking into account the energy range 4-50 keV, the sensitivity can 
be improved by 25% ( ge " siti " ity( ^°^^^ )■ 

Sensitivity ■ Eff 



a = , , V Area ■ Time (6) 

V Sensitivity ■ Eff + Background 

In this work, cr=5.5 for CRM and a=7 for ECLAIRs. Eff means the de- 
tection efficiency when the incident photons are perpendicular to the detector 



20 



Donghua Zhao et al. 



Table 1 Counting rate (counts/cm 2 /s) of each background component in various energy 
ranges. 



Instrument 


Earth 
Positions 


Background 


4-30 
kcV 


30-60 
kcV 


60-250 
kcV 


zoU-ooO 
kcV 


0.55-1 

McV 


1-5 

McV 


GRM.Nal 


Rotate 0° 


Cosmic 




0.423 


0.263 


0.018 


0.002 


0.001 


Albedo 




0.005 


0.036 


0.029 


0.013 


0.015 


Reflection 




0.001 


0.003 


0.000 


0.000 


0.000 


Rotate 180° 


Cosmic 




0.007 


0.023 


0.006 


0.001 


0.000 


Albedo 




0.108 


0.335 


0. 107 


0.036 


0.034 


Reflection 




204 


n i in 


n nm 

U.UU1 


n nnn 
u.uuu 


n nnn 
u.uuu 


GRM.CsI 


Rotate 0° 


Cosmic 






0.184 


0.026 


0.005 


0.002 


Albedo 






0.438 


0.196 


0.079 


0.095 


Reflection 






0.063 


0.002 


0.000 


0.000 


Rotate 180° 


Cosmic 




0.007 


0.023 


0.006 


0.001 


0.000 


Albedo 




0.108 


0.335 


0.107 


0.036 


0.034 


R.cflcction 




0.204 


0.119 


0.001 


0.000 


0.000 


ECLAIRs 


Rotate 0° 


Cosmic 


2.460 


0.258 


0.180 








Albedo 


0.009 


0.018 


0.109 








Reflection 


0.001 


0.002 


0.008 








Rotate 180° 


Cosmic 


0.016 


0.032 


0.108 








Albedo 


0.023 


0.073 


0.214 








Reflection 


0.119 


0.115 


0.051 









- Invalid energy range for the corresponding instrument. 



Table 2 The sensitivity of GRM_NaI t and ECLAIRs *. 





GRM.Nal 


ECLAIRs 


Earth 


FoV 


Background 


Sensitivity 


FoV 


S 


Background 


Sensitivity 


Positions 


(sr) 


(counts / cm 2 /s) 


(ph/cm 2 /s) 


(sr) 


(cm 2 ) 


(counts/cm 2 /s) 


(ph/cm 2 /s) 


Rotate 0° 


2.5 


0.41 


0.22 


2.02 


236.8 


1.10 


1.62 


Rotate 30° 


2.5 


0.42 


0.22 


2.02 


236.8 


1.11 


1.63 


Rotate 60° 


2.5 


0.43 


0.22 


2.02 


236.8 


1.15 


1.65 


Rotate 90° 


2.2 


0.44 


0.22 


1.74 


254.3 


1.06 


1.53 


Rotate 120° 


1.1 


0.44 


0.22 


0.91 


219.8 


0.74 


1.46 


Rotate 150° 


0.7 


0.53 


0.24 


0.56 


64.9 


0.57 


3.10 


Rotate 180° 


0.0 


0.61 


N/A 


0.0 


0.0 


0.64 


N/A 



' 50-300 kcV, a— 5.5, exposure time — Is. 
* 15-150 kcV, cr—7.0, exposure time — Is. 



plane in the specific energy range which is 50-300 keV for GRM and 15-150 
keV for ECLAIRs. The efficiency of ECLAIRs is computed taking the average 
GRB power-law spectrum with index of -1.64 of SWIFT as input spectrum. 
For GRM, the input spectrum is the average GRB BAND spectrum of BATSE 
(a = -1.1, /3 = -2.3, Epeak = 266 keV) [21. For ECLAIRs, the Area is S in 
Table [5] which is the weighted mean area illuminated by GRBs, considering 
the closed fraction (60%) of the coded mask, the GRB distribution in the 
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whole FoV and the Earth occultation. For GRM, the Area is 567 cm 2 which 
is the whole area of the detection crystal for two GRDs. At variance from 
ECLAIRs, GRM has no capability to locate GRBs and to distinguish the en- 
ergy deposit coordinate in the detector. The Background in the formula is in 
unit of counts/cm 2 /s. Time in the formula is equal to 1 second for both GRM 
and ECLAIRs. 

In the GRB detection-rate calculations, we assume that each Earth position 
lasts 1/7 of the whole satellite orbit period. For every Earth position, we 
calculate the GRB detection rate separately, and sum all of them together 
to get the total number of GRBs detected per year. The Earth occultation, 
the observation-time loss due to SAA passages and the trigger efficiency are 
considered. 

The GRB detection-rate calculation of ECLAIRs is based on Swift obser- 
vation. And the factor caused by the different low-energy thresholds between 
ECLAIRs and Swift is also considered. For GRM, the GRB detection-rate 
calculation is based on BATSE observation and uses the algorithm given by 
BAND [53] which was used by Dong [5], Within the assumptions, we can de- 
tecte 66 GRBs/year with ECLAIRs and 74 GRBs/year with GRM taking into 
account the crystal Nal only. Applying the same method to the crystal Csl, 
7 GRBs/year are obtained. We take the GRB detection rate coressponding to 
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Fig. 15 Evolution of the relative GRB detection rate of GRM_NaI and ECLAIRs as a 
function of the Earth positions. 



Rotate 0° as 100%. And then we get the relative GRB detection rate for each 
Earth position as shown in Figure 1151 

The main idea of the algorithm of BAND [23] is to compare the sensitivity 
to that of BATSE in the energy range 50-300 keV. It is a very good algorithm to 
estimate the GRB detection rate of GRM using the data of Nal. But for using 
the data of Csl, it is not a good choice because of its low detection efficiency 
(~5.5%) in the energy range 50-300 keV. The sensitive energy of Csl is above 
300 keV, and its detection efficiency is approximately 38% in 300-1000 keV 
taking into account the average GRB Band Spectrum of BATSE. The number 
of hard GRBs detected by SVOM can probably be increased by finding other 
methods to use the data of Csl in the energy range 300-1000 keV. This study 
will be performed in a further work. 



5 Conclusion 

In this work, three kinds of background varying with Earth positions have 
been simulated for GRM and ECLAIRs which are two high-energy instruments 
aboard the SVOM mission. In particularly, the simulations of ALBEDO with 
the model of Sazonov2007 are described in detail for its complexity and pre- 
cision. Following the background simulations, the corresponding sensitivities 
and GRB detection rates are calculated based on the given observations. 
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Concerning the impact of each background component on the instruments, 
generally, for both GRM_NaI and ECLAIRs, the cosmic gamma-rays are the 
most important background source when the Earth is outside of the FoV, es- 
pecially in the relatively low energy range. The albedo gamma-rays are always 
the most important background source in the relatively high energy range. 
The reflected cosmic gamma-rays, which are lower than cosmic gamma-rays 
in flux and lower than albedo gamma-rays in energy, have the least impact on 
the instruments, but are still an important component in the low energy range 
when the Earth enters the FoV. For GRM_CsI, albedo gamma-rays are always 
the dominant background source. 

As for the spectral models used in the simulations, the limitation lays in 
the energy range. ECLAIRs is sensitive from 4 keV to 250 keV and CRM 
is sensitive from 30 keV to 5 MeV. However, the model of Moretti2009 was 
obtained considering the data in the energy range 1.5-200 keV. The model of 
Sazonov2007 is a function of the energy which covers from 25 keV to 300 keV. 
Only the energy range of the Churazov2006 model reaches the MeV range. 
As a result, one further work on the gamma-ray background estimations is to 
find better source models for the wide energy range of SVOM. Maybe different 
models should be used for different energy ranges. 

In the calculations in Section 14.31 only gamma-ray background compo- 
nents are considered without including other kinds of backgrounds induced 
by protons, electrons, positrons, etc. as well as the background produced by 
the instrument itself. However, we cannot conclude that the value of the GRB 
detection rate is overestimated for this reason. Because there are many fac- 
tors affecting the final results in fact. For example, to compute the sensitivity 
values in Table 2 accurately, the efficiency corresponding to different incident 
directions of photons from the FoV should be considered instead of only the 
vertical direction. Furthermore, the sensitivity used to trigger GRBs on board 
will consider diverse energy ranges, integral time scales and significance lev- 
els. One assumption in our computation is that the satellite has seven states 
relative to the Earth positions and each state lasts 1/7 of the satellite or- 
bit period. This is a simplification, since the real pointing law on orbit will be 
more complex. Back to our gamma-ray background simulations, the maximum 
ALBEDO spectrum with the minimum solar modulation 0.25 GV at the high- 
est level orbit position (30°N, 287°E) was used. Accordingly, a more accurate 
estimation of the GRB detection rate needs further work and another article 
will be prepared. 
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